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Abstract 

In this work, we report our resuhs on the phase transition of a flexible homopolymer 
from a stretched chain to a compact globule. The Wang-landau method is used to study 
the thermodynamic properties of a the chain up to 512 monomers. We believe that the 
peak in the specific heat at low temperature T k, 0.05 for small chain sizes N < 100 is 
a clear evidence of the existence of metastable helical structures observed in a previous 
study. 

Keywords: Helical Structures, Wang-Landau Algorithm, Homopolymer Collapse 



1. Introduction 

In recent years, designing functional polymers at nanometers scales has gained in- 
creasing interestfl'l. However, it is challenging to control structures that have the ability 
to undergo cooperative transitions between random conformations (globules) and or- 
dered conformations (a— helices and /3— sheets). There has been a spate of theoretical 
studies dealing with this issue, commonly achieved using Monte Carlo and molecular 
dynamics simulations[2|-0]. These works have offered rewarding insights for understand- 
ing the dynamics and the phase transition of polymers and there is still a big interest 
in investigating them until today (6|. In a previous work, we have found that flexible 
homopolymers spontaneously develop helical order during the process of collapsing from 
an initially stretched conformation. We have also demonstrated that the helices are long 
lived transient states at low temperatures f?!, but with "borrowing" energy from their sur- 
rounding solvent particles, they overcome the energetic barrier and collapse into a stable 
globule (Fig. 1 shows typical conflgurations in the process of a homopolymer collapse). 
Here, we extend our study to the thermodynamic behavior of a flexible homopolymer 
using the Wang-Landau method[8]. Our main goal is to confirm the existence of the 
helical structures during the collapse, and to characterize the temperature regions corre- 
sponding to those structures. The paper is organized as follows. In Sec. II we introduce 
the model. In Sec. HI the simulation method and the parameters used to quantify the 
thermodynamic properties are presented. The simulations results are summarized in Sec. 
IV. Finally, conclusions and appropriate acknowledgments are given in Sec. V. 
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Figure 1 : Typical configurations observed during the collapse of a flexible homopolymer from an initially 
stretched chain at low temperatures, a) Initial stretched structure, b) metastable helical structure, c) 
final stable globular structure. 



2. Model 

In the simulation, we have used a simple bead-spring model. The polymer chain 
consists of N identical monomers where each monomer position varies continuously in 
three dimensions. There are two types of interactions: a harmonic interaction between 
adjacent monomers 

C/.o„aW =a(r-ro)2, (1) 
with a = lOOe/a^ and vq = 0.85cr 

and a truncated Lennard- Jones potential between non adjacent monomers 

r 4e[(a/r)i2 _ (a/r)^ + cq] for r < 2.5a 
~ \ otherwise, ^ ' 

where the constant cq is chosen such that the potential is continuous everywhere. 
Dimensionless units are used during the simulation and are defined in terms of the bead 
size a and the Lennard- Jones energy e. 

3. Method 

We have implemented the Wang-Landau algorithm to study the thermodynamic be- 
havior of the flexible homo-polymer [9j. This method is a temperature independent Monte 
Carlo technique for exploring the energy landscape. By calculating the density of states 
g{E), where E is the energy of the polymer chain, any thermodynamic observable A can 
be obtained with one single simulation for a wide range of temperatures through the 
canonical average 

Where Z is the partition function and can be determined by 

Z = Y.9iE)e-^^'^^ (4) 

E 

Using Eq. [3l one can calculate the internal energy {E)t of the homopolymer 
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and thus the specific heat C, 
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Initially, we start the simulation by taking a stretched configuration for the homopoly- 
mer chain. We set the density of states g{E) = 1.0. Then, we proceed by generating 
trial states, using appropriate "moves", and accepting them with transition probability 



where Ei is the initial energy and Ef is the energy of the trial state. 

Each Monte Carlo sweep in the simulation consists of the following moves: N diffusion 
moves, one reptation move, one single-bead crankshaft move and one pivot move. The 
moves are chosen identical to those of referencefl^. 

If the trial move is accepted, the density of states is multiplied by a modification factor 
f , i.e.g{Ef) — ^ g{Ef) x / , and a histogram H{E) is also updated, H{Ef) H{Ef) + 1. 
If the trial state is rejected, the same procedure is followed for the initial energy Ei. 
This process is repeated until the histogram H(E) is sufficiently flat, at which point the 
histogram is reset to zero, the modification factor is reduced by / ^ ^/J , and the 
random walk continues. The initial value of the modification factor is set to / = and 
is eventually reduced to a value which is /gj^^j ~ 1- The fiatness of H{E) is defined as its 
minimum divided by its average, and in this study, the flatness criterion is set to 0.6. To 
prevent overflow during the simulations, the logarithm of the density of states ln[g{E)\ 
is used instead of g{E). More details about the Wang-Landau algorithm can be found 
in reference Q 

In this study we have defined the energy range per monomer E/N = [—3.75,3.0]. 
Due to overflows in the computation, the minimum energy boundary is far from the real 
E^-nin/N = — 7.75e, found in a previous workQ. Pushing the minimum energy to the 
lower limit needs more computational efforts and has to be explored in the future. The 
bin resolution is taken dE = 0.1. All the results obtained represent averages over 20 
independent runs. 

4. Simulation results 

The temperature dependence of the specific heat for chain lengths smaller then A'' = 
100 are shown in Fig. ^ The large peak around T « 0.4 corresponds to the coil-globule 
transition. 

At lower temperatures, for N — 64, a sharp peak forms around T « 0.05, indicating 
another transition. We believe that this transition corresponds to the formation of helical 
structures. The peak becomes smaller for N = 96. 

Fig. [3] also shows the temperature dependence of the specific heat but for chain sizes 
greater than N = 100. At low temperature, for chain lengths A^ = 128 and A^ = 256, 
we observed large shoulders. This can be explained by the fact that it is difficult for 
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Figure 2: Specific heat versus temperature for cliain lengtlis A'^ = 64, 96. Two features appear: a 
large peak at temperature T ^ 0.4 representing the coil-globule transition and a sharp peak at low 
temperatures around T 0.05 corresponding to the formation of helical structures. 

helical structures to form correctly for longer cliain sizes. Those shoulders are rapidly 
vanishing for greater chain sizes {N = 512) because of the metastable nature of the 
helical structures. 

5. Conclusions 

We have studied the thermodynamic behavior of a flexible homopolymer using the 
Wang-Landau method for chain sizes up to iV = 512. Our results show two kind of 
transitions. The coil-globule transition is observed around the temperature T sa 0.4, for 
all the chain sizes studied . At low temperature T 0.05, a second transition appears 
and depends on the homopolymer chain size. For small chain sizes N < 100, a sharp 
peak appear corresponding to a long lived helical structures. However, for N > 100 we 
observe a shoulder that become small as the chain size increases. We believe that for 
infinite chain sizes, this transition vanishes and that the formation of helical structures 
occurs only for finite chain sizes. Although, the potential model used for the fiexible 
homopolymer is different from the one used by Seaton et al. [11] , our results tend to have 
some similarities. 
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Figure 3: Specific heat versus temperature for chain lengths = 128, 256 and 512. These chain lengths 
offer two features: a large peak in Cv indicative of the coil-globule transition and a shoulder indicative 
of rapidly vanishing metastable states. 
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